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Biographical data will not be published unless there 
is a demand for it. This is due to pressure for the pub- 
lication of material which seems to be of greater interest 


to our readers. 
Editor. 


EDITORIAL NOTE - 


Readers are invited to send us their selections of 
remarkable situations in Mathematics that can be explained 
at the freshman level. 

H.V. Craig, Editor. 
(Answers, continued from page 57) 
kp (A)ff - (A)f 4 = ph), f & f= ‘uayy *Ap(4),f = *P 
pue (4)f =* uaya ‘(*),_f = £49] *xp(*),_f [= I ‘oz Vv 


> £ - (yt) (G&L) st seaano 
S 


ay} UseMQ0q Bae 9yQ ‘e_Zueqoe1 Butqrs9swnd1t9 94 Fo ATeatQIeds 
-91 7/[ pue ¢/[ ere OTqno eB pue eyoqeied & sprsuT eere 943 9dUTS 
"(9‘Z) Pue (Qg‘Q) 3% 39eST9qUT SaAINd 9yq UOTIDedsUT AY *EOZ V 


(Sp = (*)f queys 
-uod 2943 eq ysnw 4t ‘spotied quapuedepuT om] sey (x)f a2UTS “goz¢ V 


*** ‘pe ‘pz ‘py ST uoTSseilZo0id ayy ‘a90us}] *p = D ALTNSeT 9yy OF Spee] 
UOTITPUOD UPATS 9yQ UT SUI9q YI(W + U) pue YQw ‘YQU ayy FO SENTRA 
aya Surynqrysqns *** ‘pz7+ vo ‘p +o ‘vy uotssaiz0id 9y44 ul ‘20z¢ ¥ 


*smoT Toy Aqtypenbout ayy ‘eoueyq “E/T = 2 
= A=X Ueya sInds0 wnwtuTW ITayy ‘snyy *auo syenba Aytypenbout 
242 JO apts puey azeaT ey UO SIOJDeF Bay. FYI Fo wNS ZY] “90g V 

*punold ay) 09 PATIBT AI SaTTw 
[Z JO Ite 943 OF SATIBTAI SeTtW GZ JO ‘ZutuINy sIOJeq ANOY y% S2TTF 
auetd aya ‘ar1ofaray] *3Sar 03 Te pue uooTTeq ayy Zutdutsg werqoid 
afoym ey 07 styy Atddy ‘ydw yt st putm ay fo peeds ay, “ggg V 
*3[Nsel paitsap ayy 393 om 


‘szoutw Ac Zutpuedxy *g = os - ° = f+yroir.ey usu 
' ee | tp 4d °O I 0 ee I 





THE RANDOM SIEVE 


by David Hawkins 


Cramer (1) discusses a simple stochastic model for the distribution 
of primes. Let there be a sequence of independent trials of an event 
S, with 


(1) Pr(S_) = I/log n 


The numbers P, , Poy see oe eee, Of which P. is the mth value of n 
for which S| occurs, will then with probability one have a limiting 
density 1/log n, like the primes. A number of other consclusions follow 
from the strong law of large numbers, for example that with proba- 
bility one 
(2) lim sup Perr fe 1 
a-© log “P, 

The model is however quite artificial in that it has the prime 
number theorem built into it ad hoc and in its assumption of inde- 
pendence. Primes are not independent in a statistical sense. The 
occurrence of an unusually long run of composite numbers leads one 
to expect a compensating increase in the number of primes later on. 

A more natural stochastic model is that of the random sieve. In 
the sieve of Eratosthenes we sieve out the multiples of every number 
which is not,a multiple of some earlier sieving number. We define 
the random sieve as follows: Check the number 2, and then with pro- 
bability % strike out each subsequent number. If P, is the first 
number not stricken out, check it and strike out each number there- 
after with probability 1/P,. P, is the next number not stricken out, 
and we use 1/P, as the probability with which to strike out each 
subsequent number, etc. The set of all possible sequences of numbers 
checked (sequences of random sieving numbers) contains the sequence 
of primes as one of its most probable members. Another typical member 
is the sequence of ‘‘lucky numbers’’ (2). S, now stands for the pro 
position that n is a sieving number, and we have the following 
recurrence relation: 


Pr(s,)? 


n 


(3) Pr(S, , ,)=Pr(S,) - 


Proof: Let T, be the contradictory of S . Let S* be the proposition 
that n is not sieved out by any sieving number less than n - 1, and 


let 7* be its contradictory. Then 
1 
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T, * tf, 1 VST R41 v Salas rast 


n+1 n°’ nt 


Pr(T.T,,,) = Pr(T,)*, Pr(S,T%,,) = Pr(S, Pr(T,), 


n+1 


Pr(S,7,,3S24,)= PriS,)?/n. 


noi 


Since the alternatives are mutually exclusive, the recurrence relation 
follows. 


We solve (3) as follows: 

Substitute g, = 1/Pr(S,), obtaining 
a Se 
n- 1/g, 


4 
(4) ur *e* 


which may be solved recursively, using g., = 5 1/n =L(n) as the 
first trial solution. This gives g, = L(n) + 0(1) and thus 


(5) Pr(s} ~ 1/log n 


Thus the random sieve gives asymptotically the result assumed 
ad hoc by Cramer. The events S are not, moreover, independent. In 
fact, it is obvious that 

(4) PERT, «os T, gig? COL 

which shows that the interdependence is negative, in conformity 
with our comment about the primes. The result does not affect the 
validity of (2), however, unless it strengthens it. For if in (4) 
we put r = c(log*n), it is easy to see that the occurrence of a 
run of that length without sieving number is asymptotically of 
probability 1/n°. From the convergence of 51/n° for c>1 it follows 
that with probability one the number of runs of such length is 
finite. Hence log*P, is almost certainly an upper bound, from 
some P. on, to the interval P,,, - P,. Because of interdependence 
it is more difficult to prove that this is false for any c > l. 

It is not difficult to define and solve various random-sieve 
problems analogous to those of multiplicative number theory, for 
example the relative frequency of pairs of sieving numbers separated 
by a given interval, the expected number of “divisors’’ of any 
number, etc. 

In the case of the lucky numbers 1, 3, 7, 8, 13, 15, ..., we 
start with odd numbers and sieve out first every third odd number, 
leaving 7 as the next lucky. Then we sieve out every seventh of 
the remaining numbers, etc. In this case nothing is known about the 
asymptotic density, so we are in the position, say, of Gauss viz 
a viz the prime distribution. A randomization of this process is, 
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however, just the random sieve again, and from this fact we im- 
mediately conjecture that the asymptotic density of luckies is 
1/logn rather than, say, B/logn where P's," 

There is one direction in which the random sieve may facilitate 
something more than conjectures. The distribution of numbers prime 
to the first m primes, or that of pairs of such numbers separated 
by a constant interval (e.g. 2) is in all likelihood more regular 
in a certain sense than the corresponding random sieve distribution, 
and if this is true it implies a number of results somewhat stronger 
than those that have been obtained, e.g. the infinity of twin 
primes. 

Suppose that the first m random sieving numbers are the first 
m primes, P. =p, r=l1, 2, ... m. Then the distribution of the 
number of numbers not sieved out from a sequence of N consecutive 
numbers > p, is given by the binomial distribution, with proba- 
bility Q = (1 - 1/2)(1 - 1/3) ... (1 - 1/p,). This same distri- 
bution may be expressed in a different way as the sum of 2" random 
variables (non-independent). Let N(i) be the number of numbers 
sieved out by the sieving number p;, N(i,j) those sieved out both 
by p; and p., etc. Those not sieved out will be given by the well- 
known combinatoric formla 


(7) N-ZN(i) + ZT Mi,j) ... +(-1)" N(1,2,3,...m). 
i i,j 

For the sieve of Eratosthenes, on the other hand, we have a 
precisely similar expression. In this case the distribution of 
numbers prime to 2.3.5. . p, =K, is periodic modulo K,, and we 
pick our sequence of N numbers, N << K,, at random from a period 
of length K,. We can calculate moments for the distribution of the 


number given by (7), both for the random and the Eratosthenes case. 
If as seems likely we can prove that moments of even order for the 
Eratosthenes distribution are smaller than the corresponding moments 
of the random sieve distribution, then it will follow that the 
longest interval between numbers prime to K, is of the order of 
Pp, log p,, a stronger result than has been obtained by other 
methods. A similar argument applied to twins prime to K, would, if 
valid, sstablish the infinity of twin primes. It is not hard to 
prove the inequality for second moments, but the problem of proof 
for moments of order 2k remains. 
*e+ 8 @ @ 

* Subsequently verified by W.E. Briggs and the author, and inde- 
pendently by P. Erdos. See *The Lucky Number Theorem" to be published 


in this magazine. 





References: 
(1) H. Cramer, Acta Arithmetica, Vol. 2, 1937, pp 23-8. 


(2) Ulam et al, On Certain sequences o futougrs defined by Sieves 
Mathematics Magazine Vol. 29 (1956) pp 117-122. 
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A Slide-Bee 
James G. Dyhikowski 


Recently there appeared in this magazine an article entitled 
“Dig That Math:’ The article was humorous and somewhat sarcastic. 
It implied that mathematics would need a radical change before 
mathematics teams would he formed. I would like to quote an article 
datelined Chicago, March 19, 1957 and entitled,“Mathematics is a 
New ‘Sport’? 

“A group of engineers who are worried about the nation’s shortage 
of engineers have devisedaslide rule competition for youngsters. 
They hope it will ‘make mathematics as popular as sports.’ 

“The first of a series of planned ‘slide-bees’ was held yesterday 
in suburban Wheaton High School. Seven-member teams from seven high 
competed, each armed with a slide rule. 

“The contest was run off like a spelling bee, with contestants 
given mathematical problems and those answering them fastest piling 
up points towards the championship. 

“Among questions posed was this one: 


How soon will the Russians overtake our supply of en- 
gineers, if they turn out 81,000 engineers every y ear 
to our 28,000 engineers? Assume that we now have 642,000 
engineers, against 396,000 engineers in Russia, and that 
deaths and retirements are disregarded. 


** The answer to this was 4 years and 7 months. 
“The winner was George Guerin of suburban Hinsdale High School. 
**Guerin’s prize: A slide rule. 


The only ambition which I have, is that this movement well spread 
over the whole country and get students interested in science and 
engineering. 











CALCULATION OF A COMPLETE SYSTEM OF TENSORS 
WITH THE AID OF SYMBOLIC MULTIPLICATION 


Lewis Bayard Robinson 


Introduction 


A tensor is a generalized covariant. Complete systems of co- 
variants are given by the solutions of complete systems of linear 
partial differential equations. 

Complete systems of tensors are given by the solutions of dif- 
ferential equations of the Riquier-Saltykov type which contain 
complete systems of linear partial differential equations as a special 
case. * 

Most of us know that the German geometers write a binary form 
thus: 


(ax)? = (ax, + a,x,)(a,x, + a,x.) = a, 5%, + 2a, 5% > + a,x, ” 


We shall use the above symbolic multiplication to help usin- 
tegrate the differential equations which arise. 

Dr. T.C. Doyle has written an important paper on the tensors of 
projective diflerential geometry.?- 

The author is extending the work of Wilczynski and his pupil 
E.B. Stouffer on the covariants of linear differential equations to 
tensors and like them he rests on the work of Lie but he uses the 
German symbolic notation as an auxiliary. 

Professor Cramlet has written a work: The Invariants of an N - 
ary Q-ic Differential Form.** But his definition of a complete 
system of tensors differs altogether from that given by the author. 
Nor does he follow the method of Lie and he does not refer to Riquier. 
He bases his work on Aronhold. It is only remotely related to the 
work of the author. 

The first paper published by the author on the above subject 
appeared in the Doklady of the Russian academy (1937). Professor 
Hlavaty of Praha wrote a review of this note and remarked that the 
author had not used a tensor notation. I have with considerable 
profit introduced a symbolic notation in this work. 


1, See Riquier, Les Systenes D’Equations aux Derivees Partielles, page 502 

2. See Grace and Young, Algebra of Invariants. 

3. Doyle, Tensor Theory of Invariants for the Projective Differential 
Geoaetry of a Curved Surface. Transactions of the American Mathematical 
Society, 1944. 

4. See, Annals of Mathematics, 1930 page 134 


5 
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Part (I) 


In a previous work the author has computed a complete system of 
semitensors associated with the system of differential equations 


Ys + Pas¥t + Pr2¥d + 911¥1 + %y2¥2 = 9 


(s) }° 


Y¥> + Pot + P22¥d + G21¥1 + 222 = 0 

In the following paper he will compute a complete system of ten- 
sors associated with the above system. 

The tensors we shall consider are of order r = 1 when there is 
only one subscript. Such a tensor is I;. I; ; is a tensor of order r#2. 

We shall now set up the system of equations defining the tensors 
when r = 2. It will then be easy for anyone to write the system 
when r = 1. 

When we were computing a complete system of semitensors the de- 
terminants of the transformation was 


The determinant of the transformation of our tensors is written 
A= Oy A2 0 
So; Soe 
0 0 
The infinitesimal transformations of A are given by 
1+¢,,t P, pt 0 
Py ,ot 1+¢,,5t 0 j=1+ P, ot + $y St + ¢'dt 
0 0 1+€'8 
The transformations of the tensors are written when r=1, thus: 
ee 3 
I, .= oh Abi I, (¢ = 3.2.3) 


When r = 2, they are written 2 


1. See Robinson, Un Systeme De Riquier et Le Calcul Tensoriel 
Bulletin de la Societe Mathematique, 1940, page 129 


2. The symbol @ indicates symbolic multiplication and I; =Tii 
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= Ay last (Ay Ao; + 47,4;,) Tio +! Ay A3i+ Aa ri 113 
si A, i494 22 +( A, isi . A, joi Tos 4 A, 43.155 


We shall compute our tensors by Lie’s method of infinitesimal 
transformations. The infinitesimal transformations are: 


: 5 9 aa 
Tay + 20, ba gbt - 21 yy pot + 2 yys'dt 


eee , 4 _)s bs 
Typ ~ Tay pot + Ty ol byy + Oy 2)dt - Ty yh24dt 


bI,0- + 215.9, ,9t =~ 21, 56,,5t + 21,56 dt 
OD y3> — Taghy2et + 143(y1 + 2h2Q)dt + 1430'S 


81,4 + 1,,(2h,, + 2 ,)0t - 1,,6,,5t + 1,,5° dt 
81,,= + 21,,(¢,, + b,,)8t 


And from the above follow the formulae denoted by (F).! 


By ol Ty )by dt + Byy Ty )by dt + OB, (Ty yp, dt + Oy (Ti) d, ydt 


+ a, (Tisderst == 21 12P, o8t + 2144 $,55¢ + 27 e'8t 


®,(1,,)%,,5t ’ O1(T, 0b, 5t ’ ®,0(T,, >, 25t : 5, (T 1), ,5t 
ME"Bt = ~ Ty,y pdt + Ly phy yot + 14 ,6,,5t 
- I, Po15t + 21,4 dt et cetera. 


We write the equations analogous to system (A) in the anterior work: 
Q: 5 (T,,) = 9 
¥,; (I,,) = 0 
(i,j,k, l=1,2) 
®,,(1,,)= - 20,, $,,(1,,) = 
®,,(7,,)= 0 ®,,(1,,) = 
,.(1,,)= + 20,, $,,(T ,2) = 


®,,(1,,)= 0 5114) = 





1. Wilczynski uses infinitesimal transformations to compute complete systems 
of covariants in his book Projective Differential Geometry. 
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oe'(T, 4) =+2I,, o¢'(Ty>) = + 2I,, 
12(T,,) = 0 12143) = - I, 
®,,(I,,) = + 2I,, ®,,(,,) = +I 
(12) = 0 ®22T 5) = 
®,,(1,,) = - 21,, 5,(T43) = 
&¢'(I,,) = +215, &-:(,,) = 
®1(T2,) = 0, ; (133) = 
1,143) = + 2053, , (I3,) = 
®5(123) = +193, %2,(I33) = 
$21(1,,) = - Igy $21,(133) = 
e'(I5,) = + To, .'(I3,) = 0. 


We derive the above from formulae (#) by equating the coef- 
ficients of the Pi; and £'. 


We know how to write the expressions, Q; (F), ¥,;(F), ®, ,( RB 


for the author has given them in the anterior work. But we have 
introduced a new symbol o,(F), which we write 
s..(F) = , OF , ¥ OF or 
a - oe a,r” @ = “Ae s 
°% °%e omy Pe 


OF , aw = 
= "is = 


where @ = 0, er?- 4J. 


Wilczynski has computed this seminvariant. See Projective Dif- 
ferential Geometry, page 110. 

We observe that of the above equation Wilczynski uses only 
the terms 


‘ bia? OF 
_ e— -5e onan 
306 30’ 


since he does not for his purposes require the other terms, when 
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he is computing only invariants. ! 
And to the equations of the above system we add one more set: 


ol ol ol ol ol 2 
ki =, Akl 2. moet. & 


‘ 


> $+ - 
* ey 7 2y, Py 2 36 


&-n (I, ) =- y 


All the equations of the above system we will denote by (A). 
Since this system (A) represents the infinitesimal transformations 
of the I,, the totality of all the solutions of (A) gives us the 
complete system of tensors. 
Now let 


Zy= 29, +%Py + YoPr2 
Z, = 2y, + YsPa1 + Y2P22 
Substituting these new variables in, (F) and d-"(F) we obtain 


5e' ms 0 
Yr ae 


36’ 
These are the equations which correspond to s (A) 
anterior work. Let us now set up the equations which correspond Lo 
system (B) in the anterior work. We write: 


Q..(F) = 0 


F 


¥,;") 


Oe (F) = 


oF OF OF 


®,.(F) + I, ,.— + 27,, — +1,.— + 21,. — + 2I.. —=-09 
1 2 22 1 2 33 
: : Fy. F 59 ‘ 981, 4 F 55 OI; 





1. See Projective Differential Geometry, page 121, formula (103). We 
replace 8, by © 
2. See the just cited formula. (103). 
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oF 


A S . 


+ + 2 


+ oF 


oF 
®,,(F) + 214 Vy + JI 


— +I1,, — + 2I,, —= 
as <a 
oF oF oF 


— - 21,,——-1,,>— =? 
a “a. ~ Mes 


o,,°) -1 


Write the new variables 


Ys 


‘ y,0' 
Y,=Z,+ 
Pe ee 
After substituting these new variables in the above equations 
we obtain 


OF or OF OF OF 
~ ~ £412 5, ~ 422 47 “23 OF 
oY, or 25, ly, 


oF oF OF oF 
—+I,,——+2I +I 


ay. 12 ar 22 ar 13 ry A 


1 12 22 13 


—— +I1,, ——: 
+ I 53 ae 23 °F, 


(I) indicates all five equations. 
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From the above five equations we have omitted,the differential 


a, OF OF ¥ 
coefficients <—T » ,» because we do not need them’ to 
Py Pig My 
compute tensors. (I) is a reduced form of (B). 





Also Y, and Y, satisfy 


Q,,(F) = 0 


¥,, (F) = 0 
be n(F) = 0 


which therefore vanish identically and appear no more in this work. 


If we can solve the five key equations (I) we can compute the 
complete system of tensors. 


Part (II). 


We shall begin by writing the five key equations for the case 
where r= 1, 

When we have solved these we can solve (II) by symbolic multi- 
plication. 

Lastly on page 14 we will write the formula which gives the 
complete system of tensors. 

The five key equations are 


= 
oe 
OF OF 


io 
°Y, eI 2 


oF oF 


——— ———— 2 
x a, 


2 


OF oF 


a, aI, 


oF oF OF OF 1 


RE SE ee ee 
") ®y, =,” "a “a, 


1. (II) denotes the five equations. 
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The fundamental solutions of (JI) are written 


4 
=e iy, 3, - y,1, } 


YI,-¥,1, 
“a f%: % _o-4 

I, =e D I, 
y, ¥, 


Re an 
@=@, =I? -4J 


as we have already stated. (See pages 8 and 9. 
Now 
(41, - 1)? = ¥f Ty? T, + 2 yy hel, + yf 1: I, 
= yf 1,° 1,- 2 yy tye Ty + yet yet, 
The above is unaltered by all transformations of 


ae! a 


I 2 4 


s, t 


1 z* 


But referring to foot note page 5, we replace I ;° I, by Ii; 
and write 
@* Ty HE)" = SMH? Tey - MATi2 +H Ty) = 
e* {yf Ip - 27, ¥2T 42 + yz 14,} 
We perceive the above is unaltered by the transformations of 


or I a | 


1! 12 22 


so it is a solution of (I). 

In a similar way we get other solutions of (I). 

In short the transformations of I,* I ak I, *I, and 
Ty Iy2 In are tsomorphic. 


verify the above statements we need only write 


am Ae 2 
oO + 2+ oOo, Fy +90ghd\: 1, + O30, +1, + 


. 2 . 
20,143,7 2 I, + 437, I, 


A2 2 2 
SSTEET| t 24 4074, ’ 24,4517 45 * Ay 22 * 24) A317 25 +43 1155 


etc. 
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It is evident that these two pairs of transformations are 

isomorphic. It follows that symbolic multiplication using the no- 

tation of Clebsch Gordan gives us the fundamental solutions of (I). 
For as these authors write 


we in our turn will write 


T,-J, =f; = TF 


In F;* F replace I ,° I, by T; 5° Then write 


F. oF. a F,., 
t J ty 


Clearly the F; ; are the solutions of (I). (Seetheorem page 34.) 


So we write 
Fy Fy= Fy = { o%(y,T, -y5I, )}? Oy 77» 1% 42 + ypTy) 
=F, = @*{y, 1, -y2l,{} % I -Y21;} 
OY VT. - CY V, + o¥y My. + ¥y% Ty} 


, oe 2 2 ae 4 V2 
{Y, I, YI} = VT. ~ 2447. + YT, 


4 . —4 = Vv > 
et{y, I, - y,1,}+ e*D I, * DY! ., - YT 55 
= -4 sl ba f ; 
= {¥,I, - YoIy}+ @7*°D I, =@7*D 1Y,I,, - YoTy, 
— — = g- 42 
=e “‘D I,- e@*D I, =e *D‘l,, 


We have now obtained all the fundamental solutions of (J ) 
symbolic multiplication. The author has proved that these solution 
are correct by direct calculation. 

We are now in a position to give the result at which we hz 
been aiming. 

The integration of system (I) and (II) is equivalent to the in- 
tegration of two systems of Riquier. ! Consequently, with the aid 
of our solutions of (I) we easily compute a complete system of 
tensors for the case r = 2. 

The general solution of (B) is written 
OCF iy Pyar Fyge Fags Fags F33i C,) 


1. See Requier, Les Systenes D’Equations Aux Derivees Partielles, 


Chapter XIII page 502. 
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where © is an arbitrary function and the C, are the fun- 


damental covariants of system (S) which Wilczynski has 
already calculated. But the system (A) given above represents the 
I 


first infinitesimal transformations oI.. of the tensors 
(See page 9%). Hence our tensors are the poluciene of (A). To ouliey 


(A) we — thus. ! 
let the ® | be six independent functions of F,, and C.. 


Write 
- =i 


Solve with respect to Fy, We get 

SP tencee F;, =¥, (C,) iw nes) ¢ 5 i 
where the ¥ are arbitrary functions. Write 
DF, ,, F 


D(I,;, 





M= 


The F;, are linear with respect to the I therefore M does not 


depend on the I. 


Now we solve (E) with respect to the I; ; and get 


= : kl 
I Sk>l NY 
Wj . i. Mi kl 


where Mi; are the minors of M 
The abme formulae give us a functionally complete system of 


tensors associated with system (S) when r= 2. 
In a similar fashion one can compute a complete system of tensors 


when r= 1. 


1. The process is described by Riquier in Chapter XIII of his book 





Baltimore, Maryland. 





SOME OPERATIONAL METHODS IN THE 
CALCULUS OF FINITE DIFFERENCES 


Joseph Talacko 


1. Introduction. 


Many formulas of the calculus of finite differences are very 
complicated. They can be simplified somewhat by introducing opera- 
tional symbols for numerical interpolation, numerical differentiation 
and integration, for summation and solutions of difference equations. 
Special forms may be derived readily by the use of operators in 
some operational schemes. Whenever we have equidistant intervals, the 
application of these operators enables one to memorize complicated 
formulas in several basic relations, especially for central or mean 
differences. 

In this paper a general operator, which possesses the distri- 
butive, commutative and associative properties is introduced as a 
generalization of calculus of symbols [1], [3], [8]. The standard 
operators of calculus of finite differences are special cases. 
Furthermore, this general operator may be expanded in infinite series 
and the inverse operator defined. The introduction of hyperbolic 
functions and their polynomial expansions are also demonstrated in 


this paper. 
2. The Calculus of Symbols. 


Let f(x) be a real, single-valued, continuous function in a 
closed interval a < x < b, possessing a continuous differential 
coefficient of order k. Let hy h,, hy, otk hi; H,, H,, ..., H, de 
real constants . 


Definition: An operator H is defined by the relation 
(2.1) Hf(x) =H, f(x +h,) + Haflx +h) +... +H fle +h). 
From (2.1) we get some particular cases: 
If h, = 0, H; = S,  ¢€2 23, 3, icc, tt Coe 
H=H,. 


A special operator H = 1, does not change the given function and a 
zero operator is defined to be H = 0. Two H operators are equal 
if, and only if, their application to any function gives identical 
results. Let 


15 
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K f(x) = Kf (x + k,) + oo. + Kf (x + ho}, 
(2.2) 
Lf(x)=Lyf(lx + 1,) +... +L fie + 1). 
The necessary and sufficient condition for K = L is that 


ree, «, =£,, k. l 2 oe a eee 


t u t t 


This operator possesses the distributive, commutative and associa- 
tive properties. It permits the formation of polynomials and com- 
bination of these polynomials. From the definition (2.1) we get 


(2.3) H{ f(x) + glx)} = Hf(x) + Hg(x) 
Furthermore: 


(2.4) (K + L)f(x) = Kf(x) + Lf (x) 


K{Lf (x) KL) f(x), 
(2.5) { [Lf(x)] = [KL] f(x) 


L(Kf(x)}] = (LK) f(x), 


and simply 
KL = LK. 


By induction, we may write, 
H(KL) = (HK)L 
(H + KUL HL + KL 
(H+ K)L = L(H+K). 


From the definition of the H-operator (2.1), we may write also an 
equation 


(2.6) H+X =K 


with an unknown operation X. 
Let us define the operation of division by 


(2.7) Hf(x) = glx). 


~% : ; : “Tr 
This operator H possesses the distributive properties, but it is 
not always commutative. 

Thus successive operations are indicated in reverse order. That is, 


(2.8) HH g(x) = glx). 


The operator HH~! gives the unique result 
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(2.9) HA~*=1. 


The H-operator may be generalized by introduction of the operator 
H., for which the limit is 


(2.10) lim H. f(x) = lim g(x,r) = g(x) 


r“Fo r“ro 
and simply 
H = limH, . 
r=To 
Similarly, we define 
(2.11) p (H) = Lim (by + b,H + DH? +... + BHT), 
il 
= lim b Ht 
lie So 
12 


the symbolical infinite series. It converges, if the H-operator is 
the zero operator. 


3. Special Operators in the Calculus of Finite Differences. 
Let tabulated values of the function f(x) be 
f(xy), flx,), f(x,), ihe 
and let the intervals be equidistant: 
" i-1 
For the operator H, we have 


(3.1) H f(x) = f(x +h) 


which is equal to the well known operator E. For every integer n, 
define an exponential operator 


H"f(x) = f(x + nh), 


For every integer m 


iGl<ste a — wD. 
m 


We have the relations: 
K® =H", K=H"/*, 


Of course, instead of a rational number n/m, we may have any real 
number » and the sequence 
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bq» Ho» coos bps eee 
if we write 
lim be = bh, 
r 


we may write 


Lim H’ f(x) = lim f(x + ph) = f(x +uh). 


r-~-@ r~-@ 


Because f(w) is continuous in the neighborhood of w = x+ uh, we 


may write simply 


lim H¢" = HY, 


r—o 


Special operators, with respect to the interval h are: 


a. The displacement operator E : 
(3.3) E f(x)= f(x +h). 
b. The forward difference operator A: 
A= #.1 
A f(x)=f(x +h) - f(x). 
c. The backward operator VY : 
V=E*A . 


where 


Vf(x) = flx +h) - f(x) = -Vf(x -h), 


d. The central-diference operator 6: 
5 = ph? . geh/2 

with hy h, 

=o +57 <5 =o). 


e. The mean-difference operator 


p= % (gh /2, prh/2) 


uf(x) =% {flx + ey + f(x - “a 
2 os 


f. The differential operator D : 


ws 
D=h *iin 


a-0 a 
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This operator gives the differential of the function f(x) with in- 
increment h. 
E* -1 ( ) - f(x) 
Dflal=h + him ip gy. I ore. 


a~0 a a~-0 a 





Recalling the definitions of the displacement and other difference 
operators, we may write recurrence formulas relating differences and 
differentials and find some simple formulas of interpolation, numerical 
differentiation, integration and summation. 

We may get other usable operators, by the introduction of hyper- 
bolic function. Let us first introduce the familiar operational 
identity. 

(3.9) E=l+A=e 


and the fundamental identities involving central differences, hyperbolic 
sine and hyperbolic cosine: 


D 


5 = ph/2_ peh/2 = -D/2 | eD/2 


am 


§ = 2 sinh (D/2) 


p= WEN? + B°h/2) = YreP/2 + @-D/2) 


(3.11) 
= cosh (D/2). 


Solving for the operator D in (3.9) we get 


@ 
r 


A 
(3.12) D = In(Q1 + A) “> (.})7°? — 
r 


rai 


This operator D may be written in the form 
6 ee 
(3.13) D=2 sinh” (5) =2in(5+"1+ 8 7/4), 


or after substitution 


(3.14) D =2 cosh~4u =“%1 + 87/4 . 


Many other relations may be found and various algebraic artificies 
for manipulating the basic identities may be discussed. We will show 
the usefulness of some operators which are similar to those of 


Lehmer [4], [5], and Michel [6]. 
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4. Some Formulas for Numerical Differentiation 


First, it is possible, taking powers of (3.12) step by step to 
get the well known formula of numerical differentiation of tabulated 
functions f(x), based on (3.9) in the form: 


2 ¥ n 
(4.1) pir) » cnr] 
r=1 ‘ 


The central diflerence formula offers more elegant results from (3.13) 
for n =2m 


: 2 

p(™) = [2 sinh7*(8/2)] 

(4.2) 
 - ] 
= §2a[]7 ~— + — (3/2)? + 
- 3 2 

As a particular case we have for m = 2, 


1 
p*= -*(1 - —8 7 + 
6 240 


ffn = Za +l, we can write 
pf2a+1) = [2 sinh 1(3/2)] 28+! 


;-2 
= p828+1/) _ ae (5/2)74 : 
(4.3) ; 


I ee. 
1 -— (8/2) +— + — (8/2)*-... J. 
‘aii Seniand We Yoicd J 


Consider the differentiation with respect to 5 of the equation (4.2) 
using the order as power exponent 
2m-1 


d ot 
(4.4) — D2" = 2m [2sinh7! (°/2)] - (14 87/4) = 
ds bh 


2m p2a°1 


hence d 
7-_—_—_*e— p2" 
2m dd 


In particular, we get for example 
p? = us? ~— 1/4 8° ¢ 77120 8" ~ ... 


5. Some Interpolation Formulas. 


Using (3.1) we may get the Newton forward interpolation formula, 
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for an interval of length h: 


Ce Pe Ria grate Cee + + Cates. 


1 >> ce. 
rai 


and similarly the Newton backward formula becomes 


@ 
(5.2) r= (1 +9) =-1+) (70% 

rel 
The Stirling interpolation formula, using relations (3.6), (3.10) 
and (3.11) may be written for n=2a 


D D 
(5.3) pmh/2 ~ cosh 2m 5 + sinh 2m =. 


The expansion of sinh aD and cosh aD is simple in terms of sinh D. 
Let 
sinhD =t sothat D=sinh™'t and 


sinh aD = sinh (m* sinh ‘ t) = f(t). 


Differentiating twice with respect to t, if f(0) = 0, and f’(0) =a 
we get the differential equation 


(5.4) (14 t7)f"(t) + tf'(t) - mf(t) = 0 


Using the Leibnitz theorem, differentiating r-times, for t = 0, we 


get a recurrence relation 


(5.5) fOT+20(0) = (m2 - 2) Ff) (0). 


Equations (5.4) and (5.5), with initial values of f(0), f'(0) enable 
us to obtain the successive coefficients in the Maclaurin series of 
sinh (mD). This gives us the Stirling formula: 

m(m2 - 1) m(m2- 1)(m2 - 32) 


E ®} —m sinh D + —J7— sinh’ + 7; sinh®D + ... 





Expressing sinh D in terms of half arguments, we obtain 


2m (4m?— 4) 


E** = 2m sinh (2/2) cosh(2/2) + sinh? (2/2) cosh(?/2) + 
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2m (4m?— 4)(4m - 4 ) 
; re = 2 sinh®(D/2) cosh /2) + ... 





4m? 
+ ] + oo (D/2) . ge 


finally , ie 

2 (m*- 1) (m*- 

baie, a tg og), Sa gs 

(5.6) 2! 3! 4! 

m(m2 — 1)(m?- 2?) 5 m*(m* — 1)(n* - 2?) - 
§ 5 


r 5! ee I 








If we let g(t) = cosh(mD) = cosh(m sinh~‘t) we get RBessel’s inter- 
polation formula from (5.4): 
f(0), f'(0) = 0 
plas th — pl 2mot)h/2 
n? m2(m2-22) m2 2-27) (n? -4?) 


= 1l4— sinh2D + —————- sinh*D + sinh®D +... 


2! 4! 6! 





D D 
E2(m-4)h — cosh (Qm- 1) 2 + sinh (2m - 1) ri 


2m(2m —- 2) 
= cosh (D/2) + ——— sinh*(D/2)cosh (D/2) 


2m(2m -2)(2m - 4)(2m + 2) 
4! 


sinh*(D/2) cosh (D/2)+... 





(2m ~1)2m(2m -2) 
(Se = 2st 18) 4 — = sinh3(D/2) 





(2m -1)(2m)(2m -2)(2m - 4) (2m4 2) 
5! 


sinh’ (D/2) + ... 





Finally 


p(ta-VDh/2 = 4 (m - %) 84 m(m - 1) 


iced 


§4 








m(m-1)(m-%) 3 m(m -1)(m@1)(m -2) 
57 + m 
3! 4! 


+ 


m(m -1)(m41)(m-2)(m-%) 55 
5! 





t+ see 
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6. Inverse Operators A} D7, y 7}, 


The meaning of the inverse operator, as the operator of inte- 
gration is best demonstrated by the equation 


1 
(6.1) f(x) = D7 'g(x) = S efx) dx +-C, 


which is the inverse operation of differentiation 
Df (x) = glx) = hf'(x). 


The indefinite integral D~! leaves an arbitrary constant of 
gration. So, as we find the operational expressions for nur 
differentiation, we may find operational formulas for nume.. 
integrations and summations. If’ the operators are taken between 
definite limits, the constant of integration vanishes, and we may 
write 


b 
[D-g(x)]” =f g(x) dx. 
a a 


To find Newton’s integration formula, we use the following expansion 
of the operator ); as a new operator @ ~ 4)". 


B B B 
OF <.39 2 BOs aD + 2H én. 0*> 
21 4! 61 


where B,, as known Bernouli numbers [4], [5]. The operator (e? -})" 


satisfies the equation 


B B B B, 
Ao @* 422 B+ O76 cd eV + 1 2 2D Ps ck. 
xl 4! 2 4! 


It is a reciprocal operator to A, It is not the most general operator, 
but all results of the operation (6.2) will be a subset of results 


of the operation Aw}, We may write 


B B, B 

(6.3) A~t=y-!, 2p -—2 p?+ —2p5 - 
21 4 61 

or 


D-*f(x) =1/h S f(x) dx 
B B B 
a fez pe C 4 3 rae . 
= A~*f(e) - r Df (x) + ro D? f (x) a D° f(x) Fe oA 


which is the known Euler-Maclaurin formula. In general we get 
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ae 1 
=A +o Ge ot 


and finally, the Newton integration formula 


and : (x) dx= A71f(x) : (x) . (x) 
; f(x) dx= fix + 5 fe) - so afte 


: 2f(x) = 3 f(x) 
+ me f(x) - 790 4 FUE) 4 <or 


Similarly, the operator heh represents the inverse operator, in 
the sense that it nullifies the other operation: 


(6.6) wd *flx) =4fle) +f - 1+ fe -2)4+ 2.0 4 


and for definite limits 


[s-1]° = 4f(b)+f(b-1)+ ... + fla+ 1) +%fla). 


The summation formulas may be developed by applying the procedure 
of differentiating with respect to an operator, finding the new 
operator, say y(x), and taking the limits a,b. 


7. Other Extensions of Calculus of Symbolic Operators, 


By the seme method as we have empolyed in sections 2,3 and 4, it 
is possible to define interpolation formulas in one exponential 
operator, by a method similar to the way Lehmer [6] defined the 
Bernoulli numbers. For example, the discussed Bessel’s formula may 
be defined as a primitive function, using hyperbolic function 


(m-4)D cosh(m - %4)D 
e = 


inh(m - %)D. 
cosh(D/2) cilieleiaciibas 





If we consider the other operator, say (E" - C)", [6], [8], we 
may develop a simple procedure for solutions of difference equations. 
If we consider the infinite series in some operators, through func- 
tional expressions, we get a reliable tool with may applications. 


sets *#+ ¢ + * & © £ 
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A FINITE SEQUENCE AND A CARD TRICK 


Ali R. Amir-Moez 


Suppose we take a number of cards, or we might choose, for ex- 
ample, C blank cards and write a number on each of them with 
repitition allowed. We pick a number a, larger than all of the 
numbers written on the cards. Holding the cards face up, we count 
from whatever number is on the top card, up to a. In this way we 
obtain a pile of cards. We put this pile down in such a way that 
the card from which we started counting is on the top of the pile, 
but not necessarily face up. In fact it is better to have it face 
down. We continue this counting until either we exhaust the original 
pack of cards, or until we get a remainder r such that r <a. With- 
out turning up the top cards we can find N, the sum of the numbers 
on those cards, knowing only C, the number of the cards, p, the 
number of the piles, r, the remainder, and a. 

Let MysMy » voey Mt be the numbers written on the top cards of 
the first, second, ..., p th. pile respectively, and N,, Np, ..., 
N. be the numbers of cards in the corresponding piles. It is easily 
observed that 
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=a-n,+1l, 
=a-n, +1, 


Adding these equalities we 


N, +N, +... +N, =pa- (ny +n, +... 


But ny +n, + se. N, and N,+N,+... 


Therefore N=pla+1)+r-C. 


Now we have the formula. Let us entertain our friends. Pick up 
a deck of 52 cards and let Jell, Q<=12, and K=13. Suppose a=15. 


Ask a friend to do the counting in the way described above in 
the general case. As an example, he may see 3 on the top and he 
counts 3, 4, ..., 15, paying no attention to the numbers on the 
cards after the 3. Let the next card on the top be 6. Then the 
counting will be 6, 7, ..., 15. Next he sees 2 on the top and he 
counts 2, 3, ..., 15. Let the next card be J = 11. He counts II, 12, 
.--, 15. Now a King comes on the top and he says 13, 14, 15. Finally 
the card on the top is 7, he counts 7, 8, ..., 13 and there are no 
more cards. Of course you are not watching him count the cards, in 
fact, you may be in a different room. He only tells you that he has 
5 piles and 7 cards are left over. Substituting in the formula for 
N we have 


N = 5(16) + 5 - 52 =35. 


So you say that the sum of the top cards is 35. Clearly the top 
cards in this example are 3, 5, 2, 11, 13, whose sum is 35. 

There is a short cut for getting N. We leave it to the reader 
to find it. 





Queens College, New York. 


(Miscellaneous Notes continued on page 41) 














MULTIPLE NUMBERS 


John A. Tierney and John Tyler 


A complex number is a number of the form x + iy where the com- 
ponents x and y are real and i = V-1. If we replace x + ty by x + ry 
where r? = -] we may regard x + ry as representing two complex numbers 
x + ly and its conjugate x - iy. It is interesting to consider the 
generalization to n components where r satisfies an equation of degree 
n. 


We define a multiple number as a number having the form = 
n 
x,r° where the x; are real andr satisfies the equation pats 


ay', a, real, with the latter equation having no mltiple roots. 
iz0 
c is a multiple number in that it represents simultaneously n + 1 


complex numbers. 
We consider the case n = 1 and take c = x + ry where r satisfies 


r? = pr +q_ having roots 


(1) oa oak (k 40.) 


We call x and y the components of c and x, + ry, and x, + ry, are 
defined to be equal if x, + yyy =%2 + Ty¥p and % + To¥y =%_ + ToY>- 
If x, + ry; =% 2 + Ty, their components are equal, otherwise the 
equation (x, -z,) + (y,; - y,)r = 0 would have two distinct roots. 

The norm N of c is the product of x +r,y and its conjugate 
2 + pry - gy’. 


The conjugate of x + ry is x + (p - r)y since their product is N. 
8 


x +My or Nc) =x? + (ry, + ry)xy + rytoy? =x 


To obtain a polar form for c, let c =x try=WVe" 
Then 
r,+r.,)0 
[e(r, J [e(r,)] = N= y2e' ae y2eP? 
and 
y= VN e7P/2 
Thus 


c =x +ry =V¥N e("P/2)8 


27 
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We call YN the modulus of c and 6 the amplitude of c. The mul- 


tiple number e(T°P/2)9 nas unit modulus and we write 


(2) e(r-p/2)6 - (9) + rs(@). 


To find the functions c and s we substitute rs and r, in (2) and 
solve for c and s to obtain 


(3) Mane. aw 
9 k 


and 


(4) s(8) = — sinh kO 


R 


From x + ry =VN [c + rs] we find that x = VN Cc; y-Ws and 


(5) ] ky 


Gg = . sinh”! = 


From (2) we also have [c(9) + rs(0)]" = c(nd) + rs(n?), a gener- 
alization of De Moivre’s theorem. For p=90, q=1 this reduces to 
[cosh 9 + rsinh @]" = cosh n@ + r sinhn@ and we can obtain coshn@ 
or sinh n@ in terms of functions of cosh 6 and sinh @ by expanding 
and equating components. 

To obtain the relationship between c(@) and s(@) we substitute 


r, and r, in (2) and multiply the results. 


1 
l = ¢* + (r, + rales + ry e* 
1 = c* + pes - qs?. 
To find the addition temtiibes for the functions c(6) and s(6) we 
substitute 9, and 4, in (2) and multiply. 


e (rop/2)% o(r=p/2)82 = 4 (r-p/2) (91492) 
[c(6,) + rs(6,)] [c(,) + rs(@,)] = clO, + 6,) + rs(O, + 8,) 
= ¢(6,)c(6,) + r{c(6,)s (6,) + c(6,)s (6,)] + (pr + q) s (6,)s(8,). 


Equating components we obtain 
(6) c(6, + 6,) = c(6,) ¢(@,) + qs(6,) s (6,) 


(7) ' = s(9,)c(6,) +s(6,)¢ (6,) + ps@ ,) s (8,) 





1957) MULTIPLE NUMBERS 


If in (6) and (7) we set 0, =-6, =-@ we obtain 
l=c (O)c (-@) + qs(O)s (-6) 
0 =s (O)c (-0) + s (-0) c (A) + ps(@)s (-0) 
from which 
-s(0) 


s(-@) = oP =-s(@) 
c*(6) + pe(@) s (@) = qs *(6) 





c (-0) = c(6) + ps(8). 


Using these results in (45) and (7) we have 


cd, - 6,) = c(6,)e (6,) + pe(8,) s (6,) ~ qs(0,) s (6,) 


s(G, -6,) = s (0 ,)c (6) - s(6,)c(6,). 


Continuing, we are able to derive formulas analagous to the double 
and half-angle formulas. If we compute the modulus of the sum of two 
multiple numbers we obtain a generalized law of cosines. We are able 
to differentiate the functions c(9) and s(9) together with their 
inverses and finally we can show that the components of a function 
of a multiple number satisfy the partial differential equation 
Uyy - qUxx - pUxy = 0 which reduces to Laplace’s eq. ir tne trig 


case and the wave eq. in the hyperbolic case. 
We illustrate one use of multiple numbers by solving the Diophan- 


tine equation 


let c= x + ry where r*=pri+q. If p= 0, q= -1, N(c) =x? + y? 
and c? =(x? - 3 xy*) +(3x4y - y?)r. 

Using c = ye"? and N= y? eP® it is easy to show that [N(c)]” = 
N(c"), n = 1,2,3,... . Hence for n=3 we have 


(x? + y2)3 = (x7 ~ 3xy7)2 + (3x2y - y3)? 
from which 
Z 
x? — 3xy? 
y= az*y - y? 


This agrees with the general solution of elementary number theory. 
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Angle of Inclination and Curvature 


David Gans 


The familiar notion of angle of inclination of a straight line, as 
presented in texts on analytic geometry and assumed in calculus books, 
seems so natural and perfectly straightforward that it comes as some- 
thing of a shock when we first learn that the notion involves a dis- 
continuity when applied to even the simplest of curves. If a straight 
line g meets the x-axis in a point P, the angle of inclination of g, 
we know, is the positive angle, less than 180°, whose initial ray 
extends from P along the x-axis in the positive direction and whose 
terminal ray extends along g in the direction in which y increases. 
The angle of inclination of a straight line perpendicular to the 
y-axis is taken either as 0° or 180°. 

To exhibit the discontinuity referred to, let us consider a very 
simple curve, say, one which is represented in some interval by a 
single-valued function y = f(x) possessing at least a first and 
second derivative. At a point where the tangent line to this curve 
is not horizontal it is clear that the angle of inclination of this 
line is defined, in accordance with the above definition, and, more- 
over, is continuous in the neighborhood of the point. The angle of 
inclination of a horizontal tangent line at a point of inflection is 
not defined, the above definition being ambiguous for horizontal 
lines. Let us, then, agree to define it to be 0° at a point ofthis 
type on each side of which f(x) is an increasing function, and to be 
180° at a point of this type on each side of which f(x) is a de- 
creasing function. This agreement is reasonable since it achieves 
continuity for the angle of inclination at these points. 

There remain the horizontal tangents at points of relative max- 
imum or minimum. If we proceed along the curve and approach a point 
of relative maximumfrom the left, we find that the angle of inclination 
tends to 0°, whereas if we approach it from the right this angle 
tends to 180°. For a point of relative minimum the angle tends to 
180° or to 0° according as the point is approached from the left or 
the right. It follows that, no matter which of the values 0°, 18 is 
taken as the angle of inclination at a point of relative maximum or 
minimum, there will be discontinuity at the point as far as the angle 


is concerned. 
on 
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To sum up, using symbols, if @ represents the variable angle of 
inclination of the tangent to the given curve, then a is a function 
of x: a= a(x), which is defined (or definable) for all values of x 
in the given interval, and is continuous for all these values except 
those corresponding to relative maxima and minima of f(x). 

The function a = a(x), of course, has no derivative at its points 
of discontinuity. A very large number of books on calculus, strangely 
enough, have overlooked this fact in their discussion of curvature 
in a plane. Differing somewhat in their definitions of curvature, 
these books all end up by computing it in the same way, namely, by 
finding the rate of change of the angle of inclination with respect 
to the length of arc. In the process they differentiate a= a(x), 
after writing it in the ambiguous forma = tan~'[f'(x)], then sub- 
stitute the result in the equat‘on da/ds = (da/dx)(dx/ds), and obtain 


(1) K = y"/(1 + y'2)3/2 


as the expression for curvature. They then proceed to use this formla 
to find the curvature at the vertex of the parabola z= 4ay, at the 
end of the minor axis of the ellipse bx? 4 a*y? = a2b?, at the point 
x = 7/2 of the curve y = sin x, and so forth, points for which their 
derivation of formla (1) is invalid because of thediscontinuity of a(x). 


Interestingly enough, the curvatures thus obtained at points of 
relative maximum or minimum are actually correct since formla (1) 
holds for such points despite the flaw in its derivation. This can 
be seen by noting that there is an alternative way of defining an 
angle of inclination which involves no ambiguity or discontinuity 
and which validates the faulty derivations without changing any of 
the formal steps. According to this alternative method, the angle 
of inclination, 7, of a straight line with slope # is the principal 
value of tan~'. Thus, the angle 7 is positive acute, negative 
acute, or zero according as@is positive, negative, or zero. Since 
T = tan-‘u is a continuous and differentiable function of u when 7 is 
restricted to principal values, and f‘(x) is assumed to be dif- 
ferentiable throughout the given interval, so will 7 = tan~*[f '(x)] 
be continuous and have a derivative, d7/dx = f"(x)/(1 +[f'(x)]?), 
throughout the interval. Substitution of this in dt/ds = (d7r/dx)(dx/ds) 
gives formula (1). 

Fortunately, discussions of curvature can be improved without 
taking the somewhat drastic step of changing the definition of an 
angle of inclination. All that is necessary isthat the angle 7 of 
the alternative definition be used as a measure of the directim of a curve 
without calling it an angle of inclination. No special name for it seems 
necessary. Incidentally, this use of Twould also represent a significant 
application of the principal values of the inverse trigonometric functims, 
such applications being quite scarce in introductory courses in calculus. 





New York University. 





NOTES ON CIRCULAR AND HYPERBOLIC FUNCTIONS 


William S. McCulley 


1. One of the aims of instruction in calculus is to prepare stu- 
dents for solving differential equations. An important class of dif- 
ferential equations consists of ordinary equations, of second and 
higher order, with constant coefficients, which describe oscil lating 
and non-oscillating electrical and mechanical systems. The former 
are quite naturally described by oscillating functions, i.e., cir- 
cular sines and cosines. The latter are also quite naturally described 
by non-oscillating functions, i.e., hyperbolic sines and cosines. 

The usual syllabus for calculus courses places much more emphasis 
on developing and using the properties of the circular functions, 
relegating the hyperbolic functions to a coverage of two orthree les- 
sons. The closely parallel analytical prpperties of the two types 
of function render it feasible and instructive to treat them to- 
gether. Further, various relations between circular and hyperbolic 
functions are interesting to explore in themselves. 


2. As an aid to remembering many of the relations between circular 
functions the following diagram, seen originally in the preface to 
Wentworth & Smith’s Trigonometry, published in the early 1920's, 


serves as a useful condensation of a large amount of information. 





Figure 1. 


It is easy to verify the following rules: 


2.1 a. Functions and co-functions appear to left and right, 
respectively. 
b. Derivatives of functions (on left) are positive. Der- 
ivatives of co-functions (on right) are negative. 
2.2 Reciprocal functions appear at opposite ends of diagonals. 
2.3 Each function equals the first function following divided by 
the second function following, in either direction. 
sin 0 
E.g., tan 0 = 
cos 0 
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Each function equals the product of the two adjoining func- 
6 6 
O csc 0, 


tions. E.g., cot &=cos 
The two functions at the top of the hexagon have finite am - 
plitudes; the other four become infinite. 
The Roman numerals indicate the pairing of functions in qu- 
adratic identities. 

‘imilarity of definitions permits a similar, but not identical, 


arrangement of hyperbolic functions, exemplifying all the above 


rules except 2.1 a. 

a 2 - oe h 

ee ae 

J \ 

sinh \— — . ? cech 

i a 1% Sis 

eet 
a ene 


Figure 2. 


We shall see in section 5 how a camparison of Figure 1 and Figure 2 
yields further information about the graphs of circular and hyper- 
bolic functions. 


3. When we are given the value of one of the circular functions and 


wish to find the values of the other five functions without consult - 
ing tables, we may use the well-known device of reference triangles. 
A simple modification enables us to use a similar procedure for 
finding the values of hyperbolic functions (1). Suppose sinh 9= 3/5. 
Draw two right triangles. On the first of these triangles assign 
values to the altitude and hypotenuse according to the definition 
of the circular sine. On the second triangle, assign the same pair 
of values, but rotated clockwise one place with respect to sides. 


Thus we have 


(1) 








Figure 3. 
On the second triangle the length of the hypotenuse is obtained by 
means of the Pythagorean theorem. This value is assigned to the un- 
numbered side of the first triangle. The values of the remaining 
five hyperbolic functions may now be obtained from the numbers on 
the first triangle by applying, the definitions of the corresponding 


circular functions. 


4. The differentiation formulas for the circular and hyperbolic 
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functions and their inverses are assumed known. (See (2), (3), (4) 
“or details.) It is pertinent to observe that teaching and learning 
these formulas by pairing or “ stacking’’ them. E.g., 


cos sin 
D,(coshu) J sinhuD,u, J =-1 (circ.); J = +1 (hyp.) 


tan Diu . ; 
D (are tanh u) * , upper function - upper sign, etc. 
+u 


5. Consider two parametrizations of the unit hyperbola x*-yel: 


x =cosh ?, y =sinh?, and 
(5.1) 


x = sec 6 > y = tan 6. 


We thus establish a correspondence between the hyperbolic angle ¢ 
and the circular angle 6 expressed by 


(5.2) cosh P= sec 6, sinh ¢ = tan 6. 


These relations can be symbolized by superimposing the hexagons of 
Figure 1 and Figure 2, from which alignment we can read off the re- 
maining four equivalences: 


(5.3) tanh = sin 6, coth ¢ = esc 9, 


sech $= cos 9, esch ¢ = cot @. 


If we solve each of these, in particular (5.2), for @ in 
terms of $, or for in terms of 6, we get 

(5.4) 6 = arc tan (sinh ¢), $ = are sinh (tan 6), etc. 

It is clear from these relations that the domain (~™,©) for ¢ cor- 
responds to the domain (-'7,%n) for 9, and that ¢ = 0 corresponds 
to 0 = 0. We see that for corresponding values on the two domains 
the functions in(5.2)and(5.3)have equal function values, where they 
are defined, or else they both become infinite. Comparison of the 
evaluated derivatives of the functions in(5.2)and(5.3)shows that the 
corresponding functions have equal slopes at corresponding points 
of their domains. 

It can be verified by differentiation that dO6/d? = sech ¢ for 
each of the relations represented by(5.4a);differentiation of each 
of the six represented by 6.4b)leads to d?/dO = sec 6. That dO/dd 

d~/d6)~ ‘may be seen from(5.2a. If we reverse these diflerentia - 
tions we obtain six equivalent integration formulas for sech ? and 
and for sec 6. 

If we wish to derive such formulas by formal integration, we may 
may proceed as follows: 
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cosh? dp 
bs hdd : i 
( ) sechddd - {= = [Ss + sinh*D + arc .tan (sinh ¢) +C, 


etc. 

It is interesting to compare the number of operations necessary to 
evaluate antiderivatives of the form(5.4b)and(5.5)with the number re - 
quired for evaluating the standard antiderivatives. In the case of 
(5.4b)the standard forma is {sec 0d = In(sec 6 + tan 6)+C, 
which requires three references to tables and one addition; (5. 4b) and 
(5.5) can be evaluated with two references, or, on a vector slide 
rule, by one setting of the indicator. 

In his work published in 1830 Gudermann used the relation(5.5) 
to define the function he named the longitude. Cayley later renamed 
it the gudermannian. Thus 


(5.6) gdd =fsech ddd = arc tan (sinh¢)+ C, 
and its inverse 

(5.7) gd-'6 =f sec 646 = arc sinh (tan 6) +C. 
From the discussion following(5.4)we find gd(0)=0, gd(t@)=+ 7, 
6. From the correspondence relations cos 0 = sech ¢, sin 6 = 
tanh , we may parametrize the unit circle x? + y? = 1 in terms of 
the hyperbolic angle ¢ by x =sech¢, y=tanh¢. Then the area of 


the portion of this circle in the first quadrant is given by the 
line integral 


(6.1) f x y sech > tanh d 
1 


A w= f | 4 


dx/db dy/dp -sech dtanh ¢ sech’? ¢ 


“ff sech ode =% gd d | = W(%4n7- 0) = 


Similarly, the length of arc is given by 


(6.2) s “f (sech’d tanh* + sech'¢) tdp 5: sech ¢dd =gd o| = 


7. The functional correspondence established in (5.2)and(5.3) gives 
us another way of integrating sec’@, hence odd powers of séc 0, as 
follows: 


free 6d6 = ff cosn sech ddd = ff cosh%ae 


7.3) 
= % (sinh ¢ cosh O +t $) + C 


= % (tan 6 sec 6 + arc sihnh (tan @)) + C. 


The same functional correspondences may also be used to establish 
hyperbolic analogs of the Wallis formulas, thus 
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$n © er) 
f sin"6 dO = f tanh" d sechd do =f sind cosh" '¢ d¢; 
on. “a he 
f cos"@ dO = f sech"* bp dp =f cosh""*!¢ dd, etc. 

0 0 0 


8. As an example of the compact expression obtainable by combining 
circular and hyperbolic functions, consider the differential equation 


(8.1) Yee + 2%, + Oy = 0, 


where a and b denote constants. The auxiliary equation is 


r-,aryb= 0, 


of which the roots are 
-a + (a? - 4b)? 
2 





rs: 


The particular solution is given by 


(8.1) 
y= e**(Acosh pt + Boiak Bt), 


4 | a? — 4b | 
where a=-%a, 6=% (U(a?-4b)) , J= =o, with the 
rule: “= 


If J=+1, use hyperbolic functions; 1f J=-l, use circular functions. 


9. Additional analogies and generalizations involving circular and 
hyperbolic functions include the fol lowing: 


a. The hyperbolic analog of de Moivre’s theorem 
(cosh 6 + j sinh ?)” = cosh n$+yJ sinh nd, where j? = 41; 
b. The generalization of circular and hyperbolic functions to 


oscillating and non-oscillating Bessel functions respective! y. 


c. The generalization of circular and hyperbolic functions to 
the elliptic functions. In particular, the Jacobian elliptic 
functions sn(u|m), cn(ulm) satisfy, for limiting values of 
the parameter m, the following equations (see [5]) 


sn(u|0) = sinu, sn(u|1) = tanh u, 


en(u|0) = cos u, en(u}1) = sech u. 
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The Derivatives of the Trigonometric Functions 


M. J. Pascual 


The following method of arriving at the derivatives of the trig- 
onometric functions might appeal to some, even if for the sole 
reason that it differs from the conventional method employing the 


sinh 
- A novel feature of this derivation is the fact that the 


derivative of the sin @ is arrived at last and then as a consequence 
we may obtain the above limit. The derivation will be for © lying in 
the first quadrant and with AG positive. Obvious modifications can 
be made for other quadrants and AO negative. 

R(x,y + 
4, (x,y + Ay) 


Q 


P(x,y) 





0 








We start with the tan 6 and cot @. By definition 
tan(@ + Ad) - tan? 
D(ten 0) = lim — = 
Ae+0 46 





y+Ay y 


x x Ay 1 Ay 


= lim 


SAP ine ae 
Ao 0 dy AG x pg49 49 


where x 7 0, 9747/2. Similarly we obtain 
1 Ay 
D(cot 0) =-— lim x — 
y2Ae-0 644 


y 
where y #0, @#0. Now to obtain the lim x —, from the 


Ae=0 Ag’ 


above diagram it is easy to see that 


area OPS < area OPR < area OAR 


39 
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% (x2 + y?) AO <%x by < % [x2 + (y + Ay)?] AO 


t= x a x?+ y?. 
~0 A@ 


+7" 


D (tan 6) = ——— for x70, 64 7n/2, 


x2 


hence 


There fore 


D (tan 6) = sec? 6 
x? + y” 
D (cot 0) = ————— y ~0, 


x2 
D (cot 6) = = esc? 8 


By using the identities relating the trigonometric functions we 
easily get 


(1) D (sec 9) = sec 9 tan 9 and D (cos 9) = - sin 9 for 04n/2 


(2) D (esc @)=-csc 8 cot 8 and OD (sin 89) = cos 6 forg70. 


To verify that the formlas for the derivatives of the sin9@ and the 


cos 9 hold for these exceptional values of 0 as well, we have by the 
quotient formula that 


6 D(sin 8) - sin 6 D( 6 
D (tea 2) cos 6 D(sin 8) - sin 8 D(cos @) for 0#n/2 


cos?0 





cos 6 D(sin 6) + sin* 6 





sec? 0 
cos*6 

solving for D(sin 6) we get D(sin 6) = cos 6 for 97 7/2. Combin- 
ing this with (2) we may say that for 0 < 6 < 7/2 

D(sin @) = cos 6 
Similarly by using D(cot 9) for @ 7 0 we obtain D(cos 6) = - sin 0 
for €7 0 which combined with (1) yields 

D(cos 6) = -sin for 0< O6< 1/2. 


" _. «ier ; _ sin(@th)-sin 6 
Finally to evaluate the lim we have D (sin 9) = lim 
h-0 OA h~0 h 





" . . sin h 4 sin h 
which for 6 = 0 gives us_ cos 0 = lim or im =], 
h-~0 h h70 A 
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ON A CHARACTERIZATION OF ORTHOGONALITY 


Waleed A. Al-Salam 


D. Dickinson [1] gave the following theorem: 
Theorem: For a set of polynomials p,(x), n = 0,1,2,... where each 


Pp, 1s of degree precisely n, to have the property 


(x) p,(-x) + po 


a - 
n Pre 1 
for some a, #0, n=1,2,..., it is necessary and sufficient that 


there exists a recurrence relation of the form 


(x) 0 (n=2,3 ) 
nPne1'*?> U, 7 Saws 


p,{x) - xB, Pn 1X) t¢ 


and that p,(0) # 0. 

He also proved the following characterization of orthogonality 
Corollary: For a set of polynomials p, (x), n= 0,1,2,... where each 
Pp, 18 of degree n, to be orthogonal, it is sufficient that thers 
exists a relation of the form 


a, = palx) pa. (-x) + pal-x) p,. s(x) 


where (-1)" a,a, < 0, for n2 2. 


In this note we obtain a more general criterion of orthogonality 
which actually reduces to Dickinson’s criterion as a special case 
and further includes the classical polynomials. We also obtain the 
orthogonality of a class of polynomials associated with the ortho- 
gonal polynomials. We first prove the following theorem: 


Theorem 1. Given two sequences of polynomials { f, (x)} and 
{g,(x)} n= 0,1,2,..., where f, (x) and g(x) are of degree nand 
n-1 respectively, fy # 0, g,(x) 7 0, g, = 0, then the necessary 
and sufficient condition that there exists a relation 


(1) fn ®) Sngr®) - fngr (*) B,{%) = a, #0 


n+1 


is that {f,} and {gy} satisfy a recurrence relation of type 


(2) u (x) = (A, x + B,) u(x) + Cu (x), (Cc, 4 a ek oe ee 


+1 nel 


Proof: If (f,(x)} and {gy ,(x)} satisfy (2) then it is easily seen 
that 


4] 
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f a (* Bay 16%) - fay 1°28, (*) = =C(fner (ee, (2) - fede. (2) 


= (-1)"C.C 4.6. GQ hyg, 
7 0 for all n =1,2,... 


Now assume that (1) holds. Then by Casorati’s theorem {f,} and 
{g,} are independent solutions of the difference equation 


u__ (x) fnoy') En. 1'*) 


nel 


u(x) f, (=) 8, (x) 








fae) Bna1'*) 


By expanding we find 


SuoitEnes ~ Sna Bae a 
oe n- 15n+ 1 n+ 15n-1 u i n aw g(x). 
,. 1 nel 





Une 


Now consider F(x) a oe a Ore If we eliminate Bnai*) 
by means of (1) we get 
fn-1%n . 6aSnsitneot S fatns18n-1 
f, 
Sasttn * Sani%uos 


fn 





F (x) = 





But a, does not involve x, and f,,, is of degree n + 1, then we get 
that F must be of the first degree. Hence 


a 
n 


ui,i(*) = (Aix +B) u(x) - ui. (x) 


4,21 


This completes the proof of the theorem. 
In particular if we have g (x) = f,(x) + (-1)"*1f (-x) then the 
expression in the L.H.S. of (1) reduces to 


(-1)" (f, (0) Ff, 4,¢-*) + f(-2) F412) 
and similarly the expression F, reduces to 


(-1)" Ff) fay (-*) - fragt Fn y 16 -*)) 
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from which we see that F (0) = 0 and hence B= 0. This shows that 
given a sequence {f,(x)} where f, is a polynomial of degree n,f 9 #0, 
then the necessary and sufficient condition that there is a relation 


fe) Ff ngy(-2) + f,(-2)f 414%) =a, 70 
is that {f ,} satisfy a difference equation of type 


(3) unig (%) = xAu(z) + Cou i(z), C7 0 


which is Dickinson’s theorem. Another way of stating this is that 
the necessary and sufficient consition for {f, (x)} and {(-1)"f (-x) } 
to be two independent solutions of (3) is that there is a relation 


fr Fnai'-*) + flO) fang) =q 70 


for some a,. 

Now the following characterization of orthogonality can be proved 
by using Favard’s theorem [2] and the previous theorems. 

Theorem 2: Given a sequence {f, (x) } where f, 1s a polynomial of 
degree n and fy # 0 then the necessary and sufficient condition 
that {f, (x) } be orthogonal is that there exists another sequence 
of polynomials {g,(x)} where g, is a polynomial of degreen - I, 
Bo = 9, gfx) A 0 such that 


f,(x)e,, y's) - fra t'* 8, (*) =a, #0 and a,aZ 0. 


It is well known that if {f (x)} satisfy the recurrence relation 
(2) then we have for each p = es... 


Fae p(t) = TOM af (x) + SMa) fF, (x) 


whe re 1") (x) and S{")(x) are polynomials in x of degrees p wid p -} 
respectively and where TO") (x) = §, T(")(x) =A.x +B, sir) 0, 
and sr) (x) =C, 70. Let us call rT!) (x) and SI) (x) for each fixe 


n the Lome! polynomials of the first and second kinds associa 
with f,(x) respectively. Then we can prove the following 


Theorem 3: Givena sequence of orthogonal polynomials {f, (x) }where 
deg f, = n, and fy / 0, then the Lommel polynomials associated with 
each fixed n are also orthogonal. 

Proof: Since the {f (x)} are orthogonal, they satisfy a difference 


equation of the second order. 
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(x) = (Ax +B) u(x) - Gy, yx), 


Unet nol 
and there exists another sequence {g,(x)} such that 


fy (*) nya'*) - Fay g(*) g, (x) =a, #0 where a,a,> 0. 


Now consider the Lommel polynomials associated with f, (x). We have 
after Palama [4]. 


(4) a, 76") (x) = Snas's) ecs'*? - faap\*) Bn -1(*) 
and 


(5) aS!) (x) =f, 98, (x) - f,(x) Ens p(*) 


Hence the expression 


[ rim) (x) sim (x) - Thm)(x) Sinh(x)], 


after substituting (4) and (5), reduces to a and by theorem 


’ ’ j 3 ne 1%n+p 
2, this theorem is established, i.e., the two sequences (ri™ (x)} 


and {s'") (x)} p = 0,1,2,3,... are orthogonal. 

The following expressions of type (1) are known for the classical 
polynomials ( see Toscano, [3] ) 
T(a +n +1) 


(n+ 1)! (a+ 1) 





LEA) ¢) FOO) — Lfa)(x) PCA (x) = 


(a) (a) (a) (a) ele 
" ae 
at (x) Rov (x) = Payts) Ry TT at 


H_ (x) G..4(x) - A, (x) G (x) = an! 
n n+: n+ n 


where L (a Xx), pl2)(x), and H (x) are the Laguerre, the ultrasherical, 
and the Nitiiie polynomials respectively, whereas F(A) (x), Ri M(x), 


G(x) are the polynomials associated with them. 


oa ee ae 
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ON NATURAL BOUNDARIES OF A GENERALIZED 
LAMBERT SERIES 


Francis Regan and Charles Rust 


1. Introduction. In 1932, Feld [1]! introduced the series 
@ } n 
F(z) in 5 od ed tad (1) 


n= a n 
ae | az 


and showed that any analytic function can be represented in this form 
plus a constant. Doyle [2] determined the regions of convergence, ex- 
pansion in power series and the inversion of power series into a 
generalized Feld series, which included (1) as a special case. 

It is the purpose of this paper to investigate the Feld series for 
natural boundaries. In establishing natural boundaries for this series, 
conditions on the sequences {a, } and {b, } of (1) are determined. Pe fore 
obtaining these conditions three theorems dealing: with natura! bound- 
aries of the Lambert series 


n 
ine tT. (2) 


n=1 1-2" 


are presented and are used in broadening the conditions for natural 


boundaries for (1). 


2. As an immediate consequence of the work of Knopp [3], three theorems 
dealing with natural boundaries of the Lambert series follow. 
Theorem 1. If {c,} is a null sequence of positive real numbers, 
then the unit circle is a natural boundary of the 
Lambert series (2). 
Since {c} is a null sequence, then 
n 


> ec 


“ 


=0 


for a denumerable set of k’s. Hence the unit circle is a natural 


boundary of (2). 
We come next to 





1. The numbers appearing in [] refer to reference given in bibliography 


at end of paper. 
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Theorem 2. If {c,} is a real sequence such that 0<c <c,< C<®, 
then the unit circle ts a natural boundary of the function represented 
by the Lambert series (2). 


‘ 


Using Knopp’s method, we let z, = e277 with (k’,k) =1, be any 


rational point on the unit circle and then from the Lambert series, we 
have 

™ “~—e a 
i in + 2 : (3) 


n=i J] ~2z"” v=l ] -z* ktn 1 -z" 


c,2 


Multiply both sides of (3) by (1 -z/z,) and take the limit as z ap- 


proaches zy radially. We examine the two limits on the right. Hence 


c3"" 
lim (1 - z/z9) => > 


ZZ 0 vj —zkv 





tke 
k 


pam S 
~1 
OS tee eee 


which, if C* = y, we see that 


lim J(1 
2~2 0 


c . 
— lim > 
y>3 


On the other hand 


(] - z/z5) > c < perenne wee (5) 


k+n ] a zn 


When n is not a multiple of k, |1 - z"|> h>O, where h=1 when k=2 
and h = sin 277/k when k > 2. Hence 


js C C 
(l-2z/z,) = <- (1-0) 50 <- 
shies Gr ok ia h 
Since (4) is unbounded and (5) bounded, we have 
n 


C Az 
lim {a - z/z.) = {- 


z~Zo 1 on z® 
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Hence the unit circle is a natural boundary of the function represented 
(2). Lastly, we have 


Theorem 3. If {c, } ts a real sequence such that 0<c,< C <@ 


and if for a denumerable set of positive integers k, the series 
@ 


> c, ,/kv is divergent, then the unit circle is natural boundary of 
val 
the function represented by the series (2). 


As before we consider the two limits in (4) and (5). From (4), we 
easily get k 
. re) CE yt . 1 Cy 

"Be i col Le eet oe y" -@, 

{-1 1 - x* k yn! v 
({4], p.177), and it follows easily that the second is bounded. When- 
ever, this is true for a denumerable set of k’s,then the theorem follows. 


3. Before discussing the general Feld series, we shall state some 
theorems which are special cases, the proofs of which are obvious. 


Theorem 4. If a= a” and b= 1, then the Feld series represents 
a function which has the circle |z| = 1/|a| as a natural boundary. 


The analog to Knopp’s extension of the Franel theorem is 
Theorem 5. If the radius of convergence of = a"b z" is greater 
than or equal to 1/ lal, and if the numbers b, are such that for a 


@ bby a. . 
certain k, the series > with 1 7 0, 1, ..., k-1, converge if 
kel py + k 
for such a ka relatively prime integer k', we set z, = (1/a)e?"*k, 


then for radial approach 
ab, 2” “ 
lim {(1 - z/z,) _ See ao, 
ZZ 1 - a™z" v=1 kv 
If for a denumerably many k > (b,,/kv) 7 0, then the circle |z| = 1/|a) 
is a natural boundary of the function represented by (1). 


Theorem 6. If the radius of convergence of 2 ab z" ts greater 
than or equal to 1/|a , and the b. are real and positive, the series 


a"b.z” 


represents a function which has the circle lz | = 1/|a| as 


Lorn? 
a natural boundary when any one of the following conditions on the 6, 


is satisfied 
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when the numbers b, forma null sequence, or 


when the numbers b, satisfy the inequalities 
0<b< b. $B<®, or 


(c) when the numbers by satisfy0 <b SB<a 
and the series 5 
(b,,/kv), 


v=1 


diverge for a denumberable set of positive integers k. 


Now let us consider the series when a, =a and 6, = 1. Only when 
|a| < 1 will be investigated for if a =1 the Feld series (1) reduces 
to the Lambert series and if la | > 1, by means of the transformation 

ath yn 


az = w, we obtain > , where het <1, when n> I. 


1- a1""y" 


az” re 
Lemma. If |a| <1, then > => for jz|.<5 (< 1. 


. =waee” Oe 





This lemma is easily established by using Weierstrass theorem on double 


series. n 


az 
Theorem 7. If la| < 1, then the series > ———, represents a 
] -az" 
function having the unit circle as a natural boundary. 
From the lemma above, we readily obtain 


kv 
a"z a 


az 
lim (1-2/z.) 5  - lim {aca s => 79 
_—~? ]-az” ne 1- 2" kv 


where z, is any rational point on the unit circle. Hence the unit 


circle 1s a natural boundary. 


Returning to the general Feld series (1), we impose two conditions 
on the numbers a,. Firstly, la, | < 1 and secondly, Lim: V2, is i. 
n~ 


Now assuming that the radius of convergence of > a"b z"is greater 


than or equal to unity, the Feld series is uniformly convergent for 


n 
a,6,z 


iz| S$ 0 < 1. Expending in a power series, we have 2 
1 - az” 


*¢ a7’. 


/d 
where A. => a” 
- d/n d 


the sequence {a,b,} may be expressed in terms of the sequence {A}. This 


By means of « generalized Moebius function [2], 


by: 


function is defined as follows. 


S(1) = 
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(u/d)-1 
dju @n/u) S(d) = 0 foru > 1. 


We can show that _— the region of convergence of > a,b.z” and for 
a,o,2 
|z| <1, meme: 2 at where a,b, = x S(n/d)A,, as follows. 
n d/n 
1 -az 


/ 
1 


d 
Multiply both sides of A, "ak M0 
1 


d 
‘bg by S(n/d) and summing over 


all disvisors of n, we obtain 
1 
d/d d,d 
z S(n/d)A, =ZS(n/d) = Pa by = = Sn/d) = Pi : ‘ 
d/n d/d 1 1 d/n d,/(n/d) “1 1 


d,d)-1 d,d)e-1 
=ZS(d) = “ o ay oy. = = S(d)a",‘ ¢) ay bys 
d/n d 1 1 “1 dyd/n 1 q @% 
1fd 
n/(dd)-1 
z S(d)aq : 


(u/d}-1 


But let n/d = u, the coefficients of each a, soa, become 2 S(d)a, 1, 


u 
which from the definition equals zero for all u = 1 and equals one 


when u = 1. But when u=-l, n/d, = 1 ord, =n. Therefore 


= S(n/d) A, =a_b. , 
d/n / d nn . 
cz 
Now, expand the Lambert series (2) in a power series, we get > —*—= 
a ) ee 
= A,z", where A, bd an . 
If we now let this A, be the coefficients of the power series we ob- 
ab.z” c 2" 
tained by expending the Feld series, we will have 5 ————— = 5 —-—— 
1 -a 2" 1-2" 


where a,b, = =S(n/d) = cy.+ We now conclude that if Je, | <C <@ 
d/n d,/d. 1 
and la, | < 1, then the Feld series and the derived Lambert series re- 


present the same function within the unit circle. Since these two series 
represent the same function within the unit circle, if we multiply both 
by the same factor and take the limit as a point is approached radially 
on the unit circle, we will obtain the same limit. Hence using Knopp’s 
theorem for the Lambert series, we easily eatablish 
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Theorem 8. If the numbers a and 6b) of the Feld series are such 


that la. <landab = & S(n/d) & c, » where for an integer k, 

n d, / 1 

all k series & (fesets,. ..! fer 1=0, 1, 0, & = 1, * ha and 
ant 

uf for such a k and a relatively prime k', we set z,= e 7h , then for 


radial approach 


a bz” 
non 


ky 

.: = 

lam \) a P ; 

i~2Z ( ee Tee ge ki 
n 


- 
lf = Fe for a denumerbale set of integral values of k, the unit 
circle will 
Feld series. 


If we make use of Theorem 2, we conclude 


be a natural boundary for the function represented by the 


Theorem 9. If the numbers a, and 6, of the Feld series are such 


that la, | <landajb, = ¥ S(n/d) Sey, where 0O<cie, <C<o@, 
| non d/n d, /d “1 n 


then the unit circle is a natural boundary for the function represented 
by the Feld series. 


From Theorem 3, we have 
Theorem 10. If the numbers a and b, of the Feld series are such 


that |a,| < 1 and ab, = 7, 5/4) P 5 fd," where 0< c, < C < wand 


Cc 1 
kv ; : 
bE, us divergent for a denumerable set of k’s, then the unit dircle 


is a natural boundary for the function represented by the Feld series. 


ses e+ *#® & & & 


BIBLIOGRAPHY 


1. Feld, J.M., The Expansion of Analytic Functions in a Generalized Lambert 
Series, Annals of Mathematics, Vol. 33 ( 1932), pp. 139-43 

2. Doyle, William, A Generalized Lambert Series and its Moebius Function, 
Annals of Mathematics, Vol 40 (1939), pp. 353-59. 

Ueber Lambertsche Reihen, Journal fur Mathematik, Vol. 242 

pp. 283-315. 


3. Knopp, Renses 
and Application of Infinite Series, London, Blackie and 


(1913 


4, <==, ==, Theory 


Sons, 1946 





St. Louis University, St. Louis, Mo. 
Xavier University, Cincinnati, Ohio. 





PROBLEMS AND QUESTIONS 


Edited by 
Robert E. Horton, Los Angeles City College 

Readers of this department are invited to submit for solution problems be- 
lieved to be new and subject matter questions that may arise in study, in 
research, or in extra-academic situations. Proposals should be soneagensed by 
solutions, when available, and by any information that will assist the editor. 
Ordinarily, problems in well-known textbooks should not be submitted. 

Solutions should be submitted on separate, signed sheets. Figures should 
be drawn in India ink and twice the size desired for reproduction. 


Send all commnications for this department to Robert E£. Horton,Los 
Angeles City College, 855 North Vermont Ave.,Los Angeles 29, California 


PROPOSALS 


313. Proposed by Sidney Kravitz, East Paterson, New Jersey. 

Show that Christmas falls on a Sunday more often. than once every 
seven years. 

314. Proposed by J.M. Gandhi, Lingraj College, Belgaum, India. 

Show that 2°? = 20.C, - C3) +4(,C, - C5) + 60Cg - Cp) + o003 
the last term being n(\C_) when n is even and -(n-1) C, whennis odd. 
315. Proposed by P.D. Thomas, Eglin Air Force Base, Florida. 

Under the restriction f(c - x)= f(x) show that 
Sola + bx) f (x) dx may be written (a + (be) /2) f§ f(x)dx. 

316. Proposed by A.K. Rajagopal, Lingraj College, Belgaum, India. 

Prove that (1 + cos nO), n an integer, has a factor (1+cos 9) 
if and only ifn is odd. 

317. Proposed by Ben K. Gold, Los Angeles City Colleze. 

Prove that (e +.x)f¢**) > (¢ --x)f¢+*) for O< x <e. 

318. Proposed by Chih-yi Wang, University of Minnesota. 
, log x 


Evaluate ite then 


319. Proposed by Barney Bissinger, Lebanon Valley College, Pennsylvania 
Cos 6 Cos 20 Cos(n-1)6 


in nO 
Show that * =cos"” 16 {1+ 
sin 


for those values of 9 for which the terms are defined. 





Tacs ¥ 


Cos @ * oar 6 Cos"* 4g 
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SOLUTIONS 
Late Solutions 


274, 286, 288, 289, 290; B. Keshava, R. Pai, Mangalore, India 
286, 287, 288; Gene B. Parrish, Durham, North Carolina 
291; A.K. Rajagopal, Lingraj College, Belgaum, India 


A Planar area 


292. [January 1957] Proposed by Eugenio Calabi and Chih-yi Wang, 
University of Minnesota. 


Find the area of the region in the real xy plane such that 


|sinh x sinh y| <1. 


I. Solution by Gene B. Parrish, Burham, North Carolina. From the 
symmetry of the area for which -l < |sinh x sinh y|<1, the area in- 
tegral may be confined to the first quadrant, with x ranging from 0 
to @ and y ranging from 0 to sinh~‘csch x = In[(e* + 1)/(e* - 1)]. 
Denoting by A the first-quadrant portion of the area, the total area 
sought is . 

4A=4 Se In [(e* + 1)/(e* - 1)] dx 


Letting (e* - 1/(e* + 1) = u, then 
: 1 
4A =-8 J Inu du/(1 -u*) = 1? 


where the last integral is given in Table 108 of Bierens de Hans’, 
Nouvelles Tables D’Integrales Definies and in Edwards’, Treatise on 


Integral Calculus, Vol. II, p. 249. 


II. Solution by Chih-yi Wang. By symmetry, it suffices to find the 
area in the first quadrant. Let A be the total area of the required 
region. By using the relation sinh7‘u = log (u + /1+u? ), we have 


@ @ 
A/4 = he sinh"! (csch x) dx = J, log (csch x+coth x) dx 


1 h © e+1 

=f igen ff leg ; dx = 7°/4 
0 sinh x 0 e* — 1 ‘ 

(Pierce 520) 

Therefore, A= 7? 
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III. Solution by Eugenio Calabi. Let x = sinh~‘(tan u cos v), 
y = sinh~!(cos u tan v). Since the Jacobian 


yen 


(u,v) 


=l1, 


the area in question is the same as that of |sinusinv| <1, 
but we know this relation holds for |u| < 7/2, lv} < 7/2. 
Hence the answer is 7%, 


An Euclidean Asymptotic Construction of e 


293. [January 1957] Proposed by Raphael T.Coffman, Richland, Washington. 


Given a line of unit length, construct geometrically a line of 
length (1 + 1/n)”, where mis an integer. 


Solution by Leon Bankoff, Los Angeles, Culifornia. Consider the 
more general problem of constructing geometrically a line of length 
(1 + a)", where n is an integer and a is the length of a constructible 
line segment. 

Extend the given segment AB = 1 so that BC = a and AC = 1 + a. 
Erect perpendiculars BB’ end CC’ to AC at B andC respectively. 
Describe arc (A,AC) cutting BB’ at D, and extend AD to cut CC’ at E£. 
Then AE = AD + AC/AB = AC*. Now describe arc (A,AE) cutting BB’ in 
F, and extend AF to meet OC’ at G. Then AG = AF - AC/AB = AE + AC=AC?. 
Repetition of this procedure leads in an obvious way to the construc- 
tion of (1 + a)" for any desired value of n. For the special case, 
(1 + lfm)", let a = 1/n. 

As e = lim (1 + 1/n) we have an asymptotic construction of e. 

n~-@ 

Also solved by Dermott A. Breault, Carnegie Institute of Tech- 
nology; David J. Cartmell, The College of Wooster, Wooster, Ohio; 
Howard Eves, University of Maine; B. Keshava R. Pat, Mangalore, India; 
Chih-yi Wang, University of Minnesota and the proposer. 


Concyclic Nine Point Centers 
294. [January 1957] Proposed by N.A. Court, University of Oklahoma. 


The nine-point centers of the four triangles formed by four con- 
cyclic points taken three at a time lie on a circle. 


I. Solution by Sister M. Stephanie, Georgian Court College, New 
Jersey. Let the circle on which the four points lie be a unit circle, 
and the four points be, in a system of complex coordinates, t,, t,, 
t, and t,. The nine-point circle of triangle t,tot, will then be 
z= s,/2 + t/2 where sy = ty + ty + t,. Its center is at s,/2 and 
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its radius in one-half the radius of the unit circle. The nine- 
point circles of the other triangles are found similarly. 

Consider the circle z = S;/2 + t/2 where Sy = ty + t2 + t3+ty- 
This circle is the locus of the centers of the four nine-point 
circles, for the substitution of a properly chosen unit vector, 
e.g-, -t,, yields each of the centers in turn. It has center at 
S,/2 and radius equal to one-half the unit circle. 

Note: The point S /2 is the point of intersection of the four 

nine-point circles since it lies on each; it is also the center of 
the equilateral hyperbola which can be drawn through the four 
given points. 
II. Solution by B. Keshava R. Pai, Mangalore, India. The radius of 
the nine-point circle of a triangle is equal to half of the cir- 
cumradius of the triangle. Therefore, the radii of the four nine- 
point circles of the four triangles are equal, since the four points 
are concyclic and so the circumcircle, is one andthe same, viz 
the circle passing through the four points. 

Now, the centre locus of a family of Rectangular hyperbolas 
passing through the vertices of a triangle is the nine-point circle 
of the triangle. Consider each of the four triangles separately. The 
four nine-point.circles are the four centre loci of the four sets 
of Rectangular hyperbolas. Now, through any four points there passes 
a Rectangular hyperbola. The Rectangular hyperbola passing through 
the four given points belongs to all the four sets and hence its 
centre lies on all the four loci. Hence, the four loci, i.e., the 
four nine-point circles pass through a point. 

Since the four nine-point circles have common radii and pass 
through a common point, the centres of them lie on a circle with 
the common point as the centre and half the circumradius as radius. 

Also solved by J.W.. Clawson, Collegeville, Pennsylvania; Huseyin 
Demir, Kandilli, Eregli, Kdz, Turkey; Howard Eves, University of 
Maine; A.K. Rajagopal, Lingraj College, Belgaum, India; Chih-yiWang, 
University of Minnesota and the proposer. 


Radix Identities 


296. January 1957 Proposed by P.A. Piza, San Juan, Puerto Rico. 
Prove that the following equalities 
385 + 439 + 547 = 367 + 475 + 529 
3852 + 4392+ 5472 = 367? + 4752 + 5297 
are true not only when the six distinct 3-digit numbers are considered 
to belong to the decimal system of numeration, but also when they 


are regarded as belonging to any system or scale of numerical 
notation with any base greater than ten. 
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Solution by James A. Painter, I.B.M. Corporation, Endicott, New 
York. A three digit number abc in any base, say rr, represents 
ar’ + br +c. Hence the problem is to show: 


(3r2 4+ 8r + 5) + (4r 27+ 3r + 9) + (Sr? + 4r + 7) = 
(3r2 + 6r + 7) + (4r74+7r +5) + (5r2 + 2r + 9) 


(1) 


(3r? +8r+ 5)2+ (4r2 + 3r + 0)? + (5r2 + 4r + 7)? = 


(2) 
(3r2+6r +7)? + (472+ Ir t+ 5)? + (Sr?+ 2r + 9)? 


for r 210. 
Collecting like powers of r gives 3 equations from equation (1) 
and 5 from equation (2). 


§+9+7=7+5++9 


(8 +3 +4)r = (64+7+ 2)r 
(3444+ 5)r2 = (3 +4 +5)r? 
25 + 81 + 49 = 49 + 25 + 81 


(56 + 54 +80)r = (84+ 70 + 36)r 

(64 +30+94+72+16 +70)r? = (36 +49+4+42+40+90) r? 
(40 + 24+ 48)r3 = (36 + 56 + 20)r? 

(25+ 16 +9)r*= (9 +16 + 25)r* 


for r 2 10. 

Obviously this set of equations hold for any r. The base of a 
positional notation system must be greater than the value of any 
symbol used in the system. Since 9 is used in equation (1) and (2), 
r> 9. That is r 2 10. 

Also solved by Haward Eves, University of Maine; Harry M. Gehman, 
University of Buffalo; Erich Michalup, Caracas, Venezuela; Wahin 
Ng, San Francisco, California; B. Keshava R. Pai, Mangalore, 
India; Chih-yi Wang, University of Minnesota and the proposer. 


Construction of an Equilateral Triangle 
297. [January 1957] Proposed by Dewey Duncan, East Los Angeles Junior 
College. 


Given a point, a line and a circle in a plane. Construct an eq- 
uilateral triangle having a vertex on each of them. Determine the 
criterion for the existence of such a triangle. 


Solution by Howard Eves, University of Maine. Designate the given 
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point, line, and circle by P, L, and C. Let L’ be obtained by ro- 
tating L about P through an angle of +60°. Suppose L’ cuts C ina 
point Q. Then the perpendicular bisector of PQ cuts L in R such that 
triangle PQRis equilateral. 

A necessary and sufficient condition for a solution is that L’ 
intersects C. Because the angle of rotation may be +40°, there may 
be as many as 4 solutions to the problem. It is easy to construct 
examples showing exactly 3, 2, 1, and no solutions. 


Also solved by J.W. Clawson, Collegeville, Pennsylvania; Huseyin 
Demir, Kandilli, Eregli, Kdz, Turkey; H.M.. Gandhi, Lingray College, 
Belgaum, India and the proposer. 


An Invariant Curve 


298. [January 1957] Proposed by Huseyin Demir, Kandilli, Bolgesi, 
Turkey. 


Let y = f(x) be a curve with the following properties 


a) f(x) = fl-x) 
b) f(x) >0 for x>0 
c) f'lx) = 0 


Determine the weight per unit length w(x) at the point (x,y) such 
that when the curve is suspended under gravity by any two points 
on it, the curve will keep its original shape. 


Solution by K.L. Cappel, Philadelphia, Pennsylvania. Assume the 
curve to be suspended at two arbitrary points A and B. Let the 
weight between A and the y intercept of the curve be W. Then at 
A, the tension in the curve can be resolved into vertical and hor- 
izontal components so that W/H = tan OG or W=H dy/dx. 

Now assume the right point of support to be moved fromA to 
A’. If the curve is to retain its shape, there must be no change 
in the forces at A. This can only be the case if H is a constant. If 
ds is the length of the segment AA’, and dW is its weight,. then the 
weight per unit length will be 


W 








dW d W i d *y/dx? 


x a a 


eee eee 
| 1 + (4752 dx 1 + (27)? 
dx dx 
which can be satisfied by any curve obeying the given conditions. 
This problem is analagous to the problem of finding the optimum 
shape of a masonry arch, when the material of the arch is the only 
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load to be supported, and it is desired to have the thrust load act 
along the neutral axis in order to eliminate bending moments. 


Also solved by the proposer 


QUICKIES 


From time to time this department will publish problems which may 
be solved by laborious methods, but which with the proper insight may 
be disposed of with dispatch. Readers are urged to submit their 
favorite problems of this type, together with the elegant solution 
and the source, if known. 


Q 204. If A, B, C and D are vectors and [ABC]is a scalar triple 
product, prove that 


[BCD] A - {CDA] B + [DAB] C - [ABC] D = 0 
[Submitted by M.S. Klamkin. ] 
Q 205. An airplane with an airspeed of 100 mph, flying into the 
wind, passes over a point just as a lighter than air balloon is re- 
leased. After a while the plane turns down-wind and overtakes the 
balloon 8 miles from the point of release and one half hour after 
it was released. Assuming constant wind and no time lost in turning, 
how far did the plane fly before turning? [Submitted by Richard K. Guy.] 
Q 206. If x, y and z are positive and if x + y + z = 1 prove that 
(1 - 1)Q)/y -1)01/z-1) 28 [Submitted by M.S. Klamkin] 


Q 207. Find the Arithmetic Progression for which the nth term+ ath 
= (n + m)th term. [ Submitted by B. Keshava R. Pai.] 


Q 208. If f(x) = f(x +1) = f(x +V2) and f(0) =V2, find f(x). 
[ Submitted by M.S..Klamkin ] 


Q 209. Find the area between y = x? and y*=32 x. 
[Submitted by John M. Howell} 


Q 210. If f(x) can be integrated in finite form, show that the in- 
verse function f~!(x) can also be integrated in finite form. 


[ Submitted by M.S. Klamkin] 
ANSWERS 
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From the Preface to: 


“THE TREE OF MATHEMATICS” 


There is a great deal being published about mathematics these 
days, and that is fine; but this book is mathematics in the sense 
that it presents the epitomes of the main branches of the subject 
beginning with high school algebra and extending far into graduate 
work. 

When writing this treatise the authors gave great attention to 
making it both meaningful and understandable, an art in which most 
of them are pastmasters.The practice was to start on ground familiar 
to everyone and construct a highway, free from road blocks, through 
the wonderful world of mathematics. 

This book is, for the most part, a response to requests for 
source material from two classes of people: those who need an ever 
increasing knowledge of mathematics in their jobs, such as en- 
gineers; and others who have gone little if any beyond arithmetic 
and either need more mathematics or just want to know “what it is 
all about.” Inherent in a satisfactory response to these requests 
is an answer to the needs of teachers and students of mathematics 
and related subjects who desire to extend their horizons by home- 
study, and to do so in a minimum of time. 

In terms of class-work, the first seven chapters and perhaps 
selections from later ones would constitute about a one semester ' 
survey course in high school, while by passing swiftly over the 
first three or four chapters the entire volume could be covered 
in a two semester survey course in college. 

But classwork is always slower than effective individual study. 
The latter could reduce the above periods to months or possibly 
weeks by study at home during the evenings. ; 

Due to its broad coverage, this treatise offers an excellent 
opportunity for “reading up” on specific topics. By looking a topic 
up in the index you will find references to the page or pages where 
it is discussed or its meaning depicted by usuage. Of course,a 
similar procedure is possible with textbooks, but it would require 
some two dozen of them to cover all the topics in this treatise. 

The Tree Of Mathematics, containing about 350 pages with 85 
cuts and pleasing format and typography, will sell for the moderate 
price of $5.50 if cash is enclosed with order or $6.00 if billing is 
required. This book will be ready for delivery in about 60 days. 
Orders will be filled in the sequence of their arrival. Address: 


The Digest Press, Pacoima, California 
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